Active elastic dimers: self-propulsion and current reversal on a featureless track 
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We present a Brownian inchworm model of a self-propelled elastic dimer in the absence of an exter- 
nal potential. Nonequilibrium noise together with a stretch-dependent damping form the propulsion 
mechanism. Our model connects three key nonequilibrium features - position-velocity correlations, 
a nonzero mean internal force, and a drift velocity. Our analytical results, including striking current 
reversals, compare very well with numerical simulations. The model unifies the propulsion mecha- 
nisms of DNA helicases, polar rods on a vibrated surface, crawling keratocytes and Myosin VI. We 
suggest experimental realizations and tests of the model. 

PACS numbers: 87.17.Jj, 05.40.-a, 87.10.+e 



Directed motion without an imposed external gradient 
in a homogeneous, isotropic environment is seen not only 
in living systems [lj but also in agitated granular matter 
0, Q • Can these apparently diverse systems be under- 
stood in a unified manner? We argue here that they can, 
and present a model which applies, suitably interpreted, 
to the movement of helicases on DNA [J], the directed 
motion of macroscopic polar rods lying on a vertically 
vibrated surface Q, the crawling of keratocytes which 
contain treadmilling actin [B[ and the walking of proces- 
sive motors Q such as Myosin VI on actin filaments 0] • 
While not losing sight of the application to particular 
organisms or devices, our focus is on the general princi- 
ples governing propulsion by rectification of an unbiased 
input active noise in a homogeneous medium. 

In all the systems mentioned above, macroscopic di- 
rected motion of the center-of-mass (CM) arises via a 
coupling to internal coordinates, as a result of two crucial 
features - an asymmetrical environment for the internal 
coordinates and external energy input. Unlike in tradi- 
tional "Brownian ratchet models" [8| of directed motion, 
the asymmetry of interest in the above systems is inter- 
nal to the motile objects, and does not lie in an external 
periodic potential. Our approach is distinct from that 
of Q where the external potential plays a central role, 
and also differs from the dynamical systems approach of 
[Io| . The present model is similar in spirit to 11 , 12| but 
simpler, and differs in several important details as seen 
below. We find an unexpected range of possible behav- 
iors, especially in the dependence of the motion on the 
details of the nonequilibrium noise. 

Our model self-propelled object is a dimer whose two 
heads are coupled by a spring, in a homogeneous, dissi- 
pative, noisy environment. The damping coefficients of 
the heads depend on the relative coordinate or strain. 
The noise on the particles is made of two parts - a ther- 
mal part whose strength is determined by a fluctuation- 
dissipation relation with the strain-dependent damping, 
and a nonequilibrium or active part, with strength inde- 



pendent of the damping, which represents the external 
energy input. 

Our results are as follows: (i) The steady-state av- 
erage of the CM velocity is in general nonzero and ex- 
hibits counter-intuitive reversals of direction as a func- 
tion of the strengths and characteristics of the drive and 
the dampings, (ii) The steady state displays two other 
key nonequilibrium features: the mean internal force as 
well as the the equal-time correlation of the relative co- 
ordinate to the CM velocity are both nonzero, (iii) Ac- 
tive noise alone will not result in propulsion, even in the 
presence of an asymmetric internal potential; the strain- 
dependent damping is an essential ingredient, (iv) The 
preceding results, obtained by perturbative analytical so- 
lution of our model Langevin equations, with the coeffi- 
cient of the stretch-dependent damping as a small param- 
eter, are confirmed in detail by direct numerical solution 
of the equations of motion of the particles. 

The heads of the dimer are two point masses m„ 
i = 1,2, with positions Xi(t) at time t and relative co- 
ordinate x = X\ — X2 connected by a spring potential 
U(x), with a minimum at x m , and acted upon by viscous 
damping and noise with a nonequilibrium component. 
Thus the Langevin equations of the particles, in the Ito 
interpretation, are 



niiVi +cti(x)vi = -diU + \/2a l (x)k B T fa + y/AiQ (1) 

where Vi are the velocities of the two heads, the overdot 
indicates a time-derivative, <9; = d/dxi, oti(x) are damp- 
ing coefficients which depend on the internal coordinate 
x, the terms containing the unit-strength, independent 
Gaussian white noise sources fa (t) and Q (t) encode ther- 
mal and nonequilibrium agitation respectively, ks and T 
are Boltzmann's constant and the thermodynamic tem- 
perature, and Ai is a measure of the external energy 
input. We consider stretch-dependent dampings with 
cti(x) > to ensure positive dissipation. For calcula- 
tional convenience and concreteness we consider smooth 
spring potentials, but the qualitative results of our model 
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hold for any confining U . 

Before solving (TT]), some general features are worth 
noting. The noise-averaged internal velocity (A) = in 
the steady state, as long as U confines x so that such 
a steady state exists. Now consider the special case 
where the ai are independent of x. By inspection of 
([T]), we see then that the individual velocities (vi) = 
even for a non-centrosymmetric U(x), despite the pres- 
ence of the nonequilibrium noises Q, and even for the 
"two-temperature" [l3| case A±/ai ^ A2/ct2- Stretch- 
dependent damping is crucial to produce drift of the CM 
coordinate in this model. 

For a stiff enough spring, the dimer will explore small 
values of x so that a, (a;) ~ 70 +^%x where, for simplicity 
and with only trivial loss of generality, we have taken the 
x-independent part of the dampings on the two heads to 
be equal. Averaging over the noise in ([1]), we find 



(71+72) 70(71-72) 
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where the prime denotes differentiation with respect to 
x. This relation connects three key quantities - mean 
drift velocity, correlation of internal coordinate and drift 
velocity, and mean internal force - each of which can be 
nonzero only away from thermal equilibrium. In particu- 
lar, we see that a nonzero mean internal force (a force 
dipole [3]) is linked to xvi correlations, and that ei- 
ther leads to drift if the damping is stretch-dependent 
(t» 7^ 0)- The above relation also elucidates the manner 
in which an internal asymmetry in x leads to macroscopic 
directed motion: the drift velocity and the internal coor- 
dinate are correlated in the presence of a nonequilibrium 
driving force. 

Before analyzing our model in detail we examine 
four examples where internal frictional asymmetry and 
nonequilibrium noise lead to directed motion. 

(i) Structural studies 15[ together with the findings 
of a recent molecular simulation 0] of the PcrA helicase 
motor on single stranded DNA are of particular interest: 
its protein domains 1A and 2A contract around ATP and 
catalyze its hydrolysis which then actively stretches them 
apart. Stretch-dependent damping as in our Eq. (fTJ) is 
encoded in the fact that in the ATP-bound (free) state 1A 
has a higher (lower) barrier to motion than 2A. The peri- 
odic potential in |4[ is centrosymmetric, and serves only 
to provide the barriers that define the mobility. Holding 
the relative coordinate out of equilibrium, in this case by 
maintaining a disequilibrium between ATP and ADP + 
Pi, results in motion in the direction of 2A. 

(ii) Polar granular rods on a vertically vibrated, hor- 
izontal plate were studied in Q]. The two ends of the 
rod have different friction, so damping depends on the 
tilt, which is the internal coordinate of interest. The 
nonequilibrium agitation being uncorrelated to this fric- 
tional asymmetry, the CM of the rod translates. 



(iii) In the crawling of cells or cell-fragments driven 
by "treadmilling actin" the relevant internal coordi- 
nate is the instantaneous degree of polymerization, av- 
eraged over all filaments. This quantity is maintained 
in a nonequilibrium steady state by the balance between 
ATP-aidcd polymerization at the leading edge, and pas- 
sive depolymerization at the trailing edge. This leads 
to more focal adhesions at the front of the cell. The 
stretched cell thus detaches primarily at the rear, result- 
ing in net translatory motion. 

(iv) In the motion of Myosin VI [7| the detachment 
of the forward head is inhibited by the extension of the 
linker connecting it to the rear head. The differential 
binding of the forward head depending on the extension 
of the linker mimics our stretch-dependent damping, and 
ATP hydrolysis provides the energy source. 

Our model thus provides a unifying understanding of 
four quite distinct self-propelled systems. 

To understand qualitatively how the dimer walks, con- 
sider the case where the nonequilibrium noise and the 
stretch-dependent damping (with a simple form interpo- 
lating smoothly between T for x > and T' < T for 
x < 0) act only on one head of the dimer, say particle 
1. Let particle 2 have a fixed damping coefficient ly- 
ing between T and V . Suppose the active noise consists 
of discrete dimer-stretching events separated by intervals 
whose mean r is much larger than the relaxation time 
Td of the dimer. Then, if noise compresses (stretches) 
the dimer, particle 1 retracts faster (slower) than parti- 
cle 2, leading to translation of the CM in the direction 
of particle 1. Suppose instead the active noise consists 
of a succession of small impulses, at intervals t < ^ 
(effectively white noise as in this paper), then the dimer 
is kicked many times before it can relax. Head 1, whose 
damping increases with stretch, will accumulate a smaller 
displacement than head 2, so that the net displacement 
will be in the direction of head 2. This latter behavior is 
in fact what our calculations and numerical studies find. 
A purely equilibrium thermal white noise, with variance 
equal to 2fceT times the stretch-dependent damping, will 
of course fail to produce any net motion, because in- 
creases in the noise amplitude compensate precisely for 
the enhanced damping. 

Analytical expressions for the inchworm speed, the av- 
erage internal force and related statistical descriptors of 
the motion can be obtained in a perturbative treatment, 
expanding the damping coefficients to leading order in x. 
The value of this approach is that it elucidates the sepa- 
rate and essential roles of the x-dependent damping and 
the nonequilibrium noises, and shows the connection of 
the mean drift speed to the mean internal force and corre- 
lations of x with the CM velocity V. We choose m, = m 
and take a harmonic internal potential U(x) = \ax 2 . 
This simplifies considerably the perturbation theory cal- 
culations that follow, without losing any of the essential 
physics. We also assume that the active noise is absent 



3 



for i < and that at t = 0, the variables x, x and V are at 
equilibrium at temperature T. A formal solution for V(t) 
and x(t) can be written, using the propagators from the 
linearized version of the equations JT]), treating all the 
nonlinearities as source terms, from which we approx- 
imate the solutions to successive orders in perturbation 
theory. Since we are interested in velocities and their cor- 
relations with position coordinates, we find it convenient 
to retain inertia, particularly in the numerical calcula- 
tion, so as to avoid ambiguities in solving our stochastic 
differential equations [161 ]. Our numerical studies are, 
however, entirely in the overdamped regime. 

Numerical simulations were performed using an Euler- 
Maruyama scheme fl7j ]. For simplicity, we take in- 
dependent damping on only one of the particles, say the 
particle 1, while the other particle has a constant value 
of the damping coefficient (72 = 0). Specifically, we 
choose ai(x) = \ [(r + T') + (r - V) tanh(x/w)] , where 
w is the width over which the strain-dependent damping 
changes over between the two extreme values of T and 
r'. In this particular form, we have 70 = (r + r')/2 and 
71 = (r — T')/2w. We scale masses, lengths and times by 
m, x m and m/70 respectively and use a constant dimen- 
sionless time step of At = 10 -2 . Averages are reported 
over n = 10 6 realizations of the noise. We choose plausi- 
ble values of a — 0.05, ksT = 0.01 based on those quoted 
in 0] for the dimeric PcrA helicase and take (T + T') = 2 
in dimensionless units. To vary 71, we fix a value of w 
and change (r — V). 

The linearized version of the equations |T]). retaining 
the nonequilibrium noise, can be solved exactly. To this 
order, the noise averaged steady state values of the CM 
velocity and the internal force vanish. However, the equal 
time correlators have a distinct contribution from the 
nonequilibrium noise. In particular, the steady state 
value of the equal time correlator between the internal 
coordinate and the CM velocity is 



(xV)o = 



{A 1 - A 2 ) 
4(7p + am) ' 



(3) 



where the subscript indicates the order in pertur- 
bation theory. Notice that (xV)q is proportional to 
AT neq = (Ai — A 2 )/jo, which has a precise interpreta- 
tion as the difference in the effective temperatures of the 
nonequilibrium noise terms of the two heads constitut- 
ing the dimer. At thermal equilibrium, by contrast, all 
equal-time correlations between x and V must vanish. 
The nonzero correlation © between x and V leads to 
nonzero steady-state averages for V and x at first order 
in perturbation theory through @. Note from and 
Q that a nonzero mean drift velocity and mean internal 
force require 7$ ^ 0, i.e., stretch-dependent friction, in 
addition to the two-temperature scenario just described. 

From equation @, we notice that since x — > — x and 
Vi — > —Vi is a symmetry of the equation, (xvi) must be 
even in the 7$ . A calculation of the equal time correlator 




FIG. 1: (Color online) CM velocity as a function of the stretch 
dependent damping coefficient 71 for various AT neq with 72 = 
and w = 4.0. AT neq was changed by fixing the value of A\ = 
1.0 and varying A^. The solid lines are the corresponding 
perturbation theory calculations from equation (|4"]l. 



(xV) to next order in perturbation theory yields 

{xV} 2 = C (AT„ eg ) + C 1 (AT„ eg ) 2 (7 1 + 72 ) 2 
+ C 2 (AT„ eg )^( 7l - 72 ) 2 
+ C 3 (AT„ eg ) 2 ( 7l 2 - 7 2 ) + C 4 A 2 ( 7l 2 - 7 2 ) (4) 

where A = A\ + A 2 and the Ct are coefficients depending 
on a, 70 and m. Notice from the above equation that if 
both 71 = 72 and AT neq — 0, (xV) vanishes and hence 
there can be no motion of the CM coordinate at this or- 
der. The slope of the V vs. 71 curve at the origin is 
proportional to the value of AT neq . Tuning the value of 
AT neq changes the slope and hence the small 71 depen- 
dence of V . However at large 71 , this effect becomes small 
and the V goes as a power of 71 . Thus there must be a 
current reversal at some 71 . Equation ([JJ indeed predicts 
current reversals as a function of the 7 j depending on the 
values of the Ai. Numerical simulation results confirm 
this (Fig. [T]) and the current reversals agree reasonably 
with perturbation theory calculations. 

With the increase in Ai, the CM velocity does not sat- 
urate in our simple model. This is because the harmonic 
spring can stretch indefinitely and thus the dimer can be 
in an infinite number of states defined by the relative co- 
ordinate x. A real molecular motor, on the other hand, 
cannot consume an indefinite number of ATP units and 
thus its velocity saturates with the increase in ATP con- 
centration. Modifying U(x) so that the effective range of 
x is limited should give rise to a saturation of the CM 
velocity with increase in the input energy. 

In the presence of an external load we find a linear 
relation, to lowest order in perturbation theory, between 
the steady state CM velocity and the applied force. The 
reason for a linear relation is that the applied external 
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force does not alter the mechanism of energy uptake in 
our simple model. 

We also find a generalized efficiency 18| 77 ~ 
2m r )o{V) 2 / A to leading order in perturbation theory, in 
the absence of an external load. With 71 = 0.1, 72 = 
and Ai = 1.0, A 2 = 0.0, we find 77 w 0.2%. 

The results above are for the case of a harmonic spring 
potential and a non-centrosymmetric (Xi(x). A non- 
centrosymmetric U{x) and a centrosymmetric oti(x) also 
yields directed motion. Crucially, even in this case, the 
damping must be x-dependent for a nonzero drift veloc- 
ity. A non-centrosymetric U(x) alone does not lead to 
nonzero (V). We have confirmed this with an explicit 
numerical simulation. 



A few words comparing our model to that of Al- 
though the idea of stretch-dependent damping is present 
in these papers, our model is much simpler in the way the 
separation of equilibrium and nonequilibrium forces are 
presented. In 11[, stretch dependent damping is induced 
by activity whereas in our model it is present even in 
the absence of the nonequilibrium driving noise. It does 
not lead to directed motion because of the the strength of 
the equilibrium noise exactly balances the dissipation ac- 
cording to the fluctuation-dissipation theorem. Directed 
motion is induced by the acive noise whose strength is 
independent of the ai(x). Also in our model, equation 
([2]) explicitly clarifies the manner in which asymmetry 
in an internal degree of freedom is coupled to a macro- 
scopic coordinate in the presence of nonequilibrium noise, 
and leads to directed motion. That this is also propor- 
tional to a nonzero average internal force highlights the 
nonequilibrium nature of the phenomenon. 

A likely realization of this model is [13] in the form 
of a colloidal bead with a polymer tail. The damping 
on the CM of the polymer will depend significantly on 
its stretch, while that on the bead will not. Subjecting 
this composite colloid to nonequilibrium noise (chemical 
reactions, catalysis at its surface [2£j, fluctuating laser 
interference patterns) should cause it to drift in the di- 
rection of its instantaneous orientation. One-dimensional 
versions could be constructed using optical tweezers in a 
line trap geometry. This system could serve to test our 
theory. 

Tests should focus on ((2]) as well as the phenomenon of 
current reversal (Fig. [I]). In our model the dampings and 
their x-dependences are properties intrinsic to the dimer 
in the absence of active noise. Changing the strength 
of the active noise should change the numerators in 
leaving the denominators unchanged, thus allowing a test 
of the two equalities in ((2|). Current reversals (Fig. [T]) 
are best probed by altering the active noise levels on each 
head of the dimer. 

Several natural generalizations of our model suggest 
themselves, and will be investigated in the near future. 
These include Poisson or other active noises, dimers mov- 
ing in more than one dimension, coupled arrays or an 



elastic continuum of active particles [2l[, and coupling 
to hydrodynamic flow. 

We thank P. Perlekar, A. M. Jayannavar, D. Lacoste 
and K. Schulten for fruitful discussions. The Centre for 
Condensed Matter Theory is supported by the DST, In- 
dia. SR and MR acknowledge support from IFCPAR 
grant 3504-2. SR acknowledges a J. C. Bose Fellow- 
ship of the DST, India. This research was supported 
in part by the National Science Foundation under Grant 
No. PHY05-51164, and by the ICMR, UCSB. 



* Electronic address: vijayk@physics.iisc.ernet.in 

' Electronic address: sriram@physics.iisc.ernet.in 
at CMTU, JNCASR, Bangalore 560064, India. 

* Electronic address 



[10 
[11 



[12. 
[13; 

[14 



[15 
[16 



[17 

[is: 

[19 



Also 



madan@rri.res.in; Also at NCBS 
(TIFR), Bangalore 560065, India. 

D. Bray, Cell Movements: from molecules to motility 
(Garland Publishing, New York, 2001). 
D. Yamada, T. Hondou, and M. Sano, Phys. Rev. E. 67, 
040301 (R) (2003). 

S. Dorbolo, D. Volfson, L. Tsimring, and A. Kudrolli, 
Phys. Rev. Lett. 95, 044101 (2005). 

J. Yin, T. Ha, and K. Schulten, Biophys. J. 91, 2097 
(2006). 

A. Verkhovsky, T. Svitkina, and G. Borisy, Curr.Biol. 9, 
11 (1999). 

R. D. Vale and R. A. Milligan, Science 288, 88 (2000). 
D. Altman, H. L. Sweeney, and J. A. Spudich, Cell 116, 
737 (2004). 

F. Jiilicher, A. Ajdari, and J. Prost, Rev. Mod. Phys. 69, 
1269 (1997); P. Reimann, Phys. Rep. 361, 57 (2002). 
A. Ajdari, J. Phys. 4, 1577 (1994); I. Derenyi and T. 
Vicsek, Proc. Natl. Acad. Sci. USA 93, 6775 (1996); M. 
Porto, M. Urbakh, and J. Klafter, Phys. Rev. Lett. 84, 
6058 (2000); S. Cilia, F. Falo, and L.M. Fiona, Phys. 
Rev. E. 63, 031110 (2001); D. Dan, A. M. Jayannavar, 
and G. I. Menon, Physica A. 318, 40 (2003). 
S. Denisov, Phys. Lett. A. 296, 197 (2002). 
A. Mogilner, M. Mangel, and R. J. Baskin, Phys. Lett. 

A. 237, 297 (1998); H. C. Fogedby, R. Metzler, and A. 
Svane, Phys. Rev. E. 70, 021905 (2004). 

B. Norden, Y. Zolotaryuk, P. L. Christiansen, and A.V. 
Zolataryuk, Phys. Rev. E. 65, 011110 (2001). 

B. Schmittmann and R. K. P. Zia, Phase Transitions and 
Critical Phenomena Vol 17 (Academic Press, New York, 
1995). 

R. A. Simha and S. Ramaswamy, Phys. Rev. Lett. 89, 
058101 (2002); Y. Hatwalne, S. Ramaswamy, M. Rao, 
and R. A. Simha, Phys. Rev. Lett. 92, 118101 (2004). 
S. Velankar et ai, Cell. 97, 75 (1999). 
See A. W. C. Lau and T. C. Lubensky, Phys. Rev. E. 
76, 011123 (2007) and references therein for how to treat 
position-dependent noise without inertia. 
P. E. Kloden, E. Platen, and H. Schurz, Numerical so- 
lution of SDE through Computer experiments (Springer- 
Verlag, Berlin-Heidelberg, 1994). 

I. Derenyi, M. Bier, and RD. Astumian, Phys. Rev. Lett 
83, 903 (1999). 

R. Golestanian, Private communication. 



■5 



[20] G. Riickner and R. Kapral, Phys. Rev. Lett 98, 150603 J.-F. Joanny, F. Jiilicher and J. Prost, Phys. Rev. Lett 

(2007). 90, 168102 (2003). 

[21] This would be a realisation of a motile adhesive gel as in 



